SPECTRAL ASYMPTOTICS FOR ROBIN PROBLEMS 
WITH A DISCONTINUOUS COEFFICIENT 



c3 



Gerd Grubb 

Department of Mathematical Sciences, Copenhagen University, Universitetsparken 
| 5, DK-2100 Copenhagen, Denmark. E-mail grubb@math.ku.dk 

Ch ! Dedicated to the memory of M. Sh. Birman (1928-2009) 

a 

Abstract. The spectral behavior of the difference between the resolvents of two realizations 
A\ and A2 of a second-order strongly elliptic symmetric differential operator A, defined by 
different Robin conditions vu = bij^u and vu = &27o u i can m the case where all coefficients 
P I ■ are C°° be determined by use of a general result by the author in f984 on singular Green 

operators. We here treat the problem for nonsmooth &j. Using a Krein resolvent formula, 
we show that if b\ and £>2 are in Leo, the s- numbers Sj of (A\ — A) -1 — (A2 — A) -1 satisfy 
Sjj 3 '^ n ~ 1 ' < C for all j; this improves a recent result for A = —A by Behrndt et al., that 
yZ j 8? < 00 for p > (n — 1) /3. A sharper estimate is obtained when b\ and &2 are in C e for 
some e > 0, with jumps at a smooth hypersurface, namely that Sjj 3 /( n ~ — > c for j — > 00, 
with a constant c defined from the principal symbol of A and 62 — b\. 

As an auxiliary result we show that the usual principal spectral asymptotic estimate for 
pseudodifferential operators of negative order on a closed manifold extends to products of 
■ pseudodifferential operators interspersed with piecewise continuous functions. 

o 

O I Introduction. 

Consider a second-order strongly elliptic symmetric operator 

^ ■ (0.1) A = - dj(a jk d k u) + a u 

j,k=i 

on a bounded smooth domain O C M. n , and denote by A 1 , A v , resp. A, the realizations in 
^2(0) defined by the Dirichlet condition 70W = 0, the Neumann condition vu = 0, resp. a 
Robin condition vu — b^ou = with b real. Here 70 it = w|an, and 

n 

(0.2) vu= rij'yo(ajkdku), 

3,k=l 
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the conormal derivative, with n = (ni, . . . , n n ) denoting the interior normal to <9fi. It 
is a classical result of Birman [B62], shown also for exterior domains, that the difference 
between the resolvents of the Robin realization and the Dirichlet realization is compact 
and has the spectral behavior, for large negative A, 

(0.3) Sj((A — A) -1 - (A, - X)- 1 )] 2 /^-^ < C for all j; 

here Sj(T) denotes the j-th eigenvalue of (T*T)a (the j-th s-number or singular value of 
T), counted with multiplicities. This was shown assuming merely that b G L OQ (dfl). 

For the situation where all coefficients are C°°, the estimate was later improved to an 
asymptotic estimate 

(0.4) s J ((l-A)- 1 -(A 7 -A)- 1 )j' 2 / (n - 1 ) ^cfor j ->■ oo; 

this follows from Grubb [G74], Sect. 8 (with generalizations to higher-order operators), 
and Birman and Solomiak [BS80] (including exterior domains). 

The paper [G84] gave tools to extend (0.4) to nonselfadjoint situations (also for exterior 
domains by a cutoff technique), by showing that for any singular Green operator G on Q 
of order — t < and class 0, 

(0.5) Sj{G)j t ^ n ~ V) -> c(g°) for j ->• oo; 

here G belongs to the calculus of pseudo differential boundary operators, introduced by 
Boutet de Monvel [B71] and further developed in [G84], [G96]; c(g ) is a constant derived 
from the principal symbol g°. In fact, the resolvent difference in (0.4) is a singular Green 
operator of order —2 and class 0, when all coefficients are C°°. 

Considering another resolvent difference, J. Behrndt, M. Langer, I. Lobanov, V. Lo- 
toreichik and I. Popov showed in a recent paper [BLLLP10], on the basis of a theory of 
quasi-boundary triples by J. Behrndt and M. Langer [BL07], that when A = —A (hence 
vu = 7itt = ^2- rij'jodju) and b is a real function in Loo(<9fi), the difference between the 

resolvent of A and the resolvent of the Neumann realization A v satisfies an estimate with 
2 replaced by 3, for A in the intersection of resolvent sets g(A) fl q{A v ): 

(0.6) (1- A)" 1 - (A v - A)" 1 G C p for p > 3/(w - 1); 

here C p denotes the space of compact operators T with singular value sequences [sj{T))j^ G 
£ p ; the Schatten class of order p. (Besides real 6, also functions with a fixed sign on Im6 
were treated.) 

In the case where b G C°°(<9fi), this follows from (0.5) since (A — A) -1 — (A u — A) -1 is 
a singular Green operator of order —3, leading to the stronger estimate: 

(0.7) Sj((A- A) -1 - (A v - \)-i)f/(n-i) ^ c f or j ^ oo ; 

this was noted also in [GlOa], Cor. 8.4 and Ex. 8.5. 

The result of [BLLLP10] is more general by treating nonsmooth b. Their main result 
Theorem 3.5 is proved in a formulation where the boundary condition is "foil = G71W for 
an operator G, but it is explained in their Remark 3.7 how this can be made to include 
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general conditions 71 it = b^QU by use of the relations point of view of [BL07]; more on this 
in [BL10]. 

The main purpose of the present paper is to show spectral asymptotics estimates as in 
(0.7) for nonsmooth b. First we show (in Theorem 2.2) that an upper bound 

(0.8) Sj{(A- A) -1 - (A v - \)-i)f/(n-i) < C for all j, 

holds for any complex b G L oo (d0,); this implies (0.6), for general A as in (0.1). As a 
corollary, a similar estimate holds for (Ai — A) -1 — (A 2 — A) -1 , when the A^ are defined 
by boundary conditions vu = b^ou with bi G L 00 (dfl). 

Next, we show (in Theorems 2.5 and 3.4) that asymptotic estimates hold when the bi 
are piecewise slightly better that continuous. Since asymptotic estimates are not additive, 
we must aim directly for [A\ — A) -1 — (A 2 — A) -1 . 

Theorem 0.1. Assume that b\, b 2 and b 2 — bi are piecewise C s on <9fi for some e > 0, 
having jumps at C°° hyper surf aces. Then 

(0.9) 8j ((li - A)" 1 - (A 2 - A)- 1 ))//^- 1 ) c for j 00, 

where c is defined from the principal symbol of A and b 2 —b\. 

It suffices in fact that 61, b 2 and b 2 — bi are piecewise in H^{dQ) for some r > and 
p > (n — l)/r, see the details below. 

For the proof of (0.8) the method is, as in [BLLLP10], an application of functional 
analysis, building on a general theory of extensions (here Grubb [G68]) together with 
known facts on elliptic boundary value problems. The proof of (0.9) in the nonsmooth 
situations draws on methods and results for pseudo differential boundary operators in [G84] 
and a result on restricted kernels of pseudo differential operators by Laptev [L77, L81]. 

As an auxiliary result of independent interest we show (Theorem 3.3) that a product of 
classical pseudodifferential operators of negative order on a closed manifold, interspersed 
with piecewise continuous functions having jumps at a smooth hypersurface, has a principal 
spectral asymptotics estimate as in the smooth case. Moreover, we extend (0.4) to b G 
Loo(dfi) (Theorem 2.4). 

Some spectral estimates for resolvent differences in interior and exterior domains have 
been described recently by Malamud in [M10], and spectral asymptotics have been shown 
in [G10]; both papers treat higher-order operators but do not aim for the special bounds 
obtained here. Let us also mention that we do not here address the question of nonsmooth 
domains, as e.g. in Gesztesy and Mitrea [GM09, GM09a, GM10] and Abels, Grubb and 
Wood [AGW10], [G08], and their references. 

To keep the paper short, some introductory material found in other sources will not be 
repeated here. 

The main details of the extension theory [G68]-[G74] have been recalled and explained 
in several recent papers [BGW09], [G08], [G10a]; resulting Krein-type resolvent formulas 
are shown in [BGW09]. 

Sobolev spaces are recalled in numerous places. The basic facts we shall need on these 
and other function spaces such as Besov and Bessel-potential spaces, are recalled e.g. in 
[AGW10], Sect. 2. 

The calculus of pseudodifferential boundary operators is explained in Boutet de Monvel 
[B71] and in [G84], [G96], [G09]. 
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1. The Robin realization. 

Let O be a bounded smooth subset of R n with boundary <9fi = E, and let 

n 

(1.1) a(u,v) = ^2 (a jk d k u,djv) + (a u,v), 

be a sesquilinear form with coefficients in C°°(Q) such that the associated second-order 
operator (0.1) is formally selfadjoint and strongly elliptic. We assume moreover that a(a, u) 
is real for u G i7 1 (0) and (with c > 0, k > 0) 

(1.2) a(a, u) > c||u||f — fc||u||o, for u G H\n). 

This holds if the matrix (ajk(x))™ k=1 is real, symmetric and positive definite and ao(x) is 
real, at each 

Let b G L 00(E), and define the sesquilinear form a& by 

(1.3) a b (u,v) = a(u,v) + (&7o«, 7o^)l 2 (e)- 

Since ||7o^||| 2 ( E ) — c 'IMIi — £ \\ u \\i + C(e) \\u\\q for any e, we infer from (1.2) that 

(1.4) Reab(u,u) > ci\\u\\1 — ki\\u\\l, for u G -ff 1 (0), 

where ci < c is close to c and fci > k is a large constant. 

The sesquilinear form aj on 7 = i7 x (0) in i7 = L 2 (0) defines a realization A of A by 
Lions' version of the Lax-Milgram lemma (as recalled e.g. in [G09], Ch. 12), with domain 

(1.5) D(A) = {«£ H\Q) n D(A mSLX ) \ (Au, v) = a b (u, v) for all v G H\tt)}. 

The operator A is closed, densely defined with spectrum in a sectorial region in 
{Re A > —ki}, and its adjoint A* is the analogous operator defined from 

(1.6) al(u,v) = a(v,u) + (bjou, 7o^)l 2 (e)- 

In particular, when b is real, A is selfadjoint. 
It will be useful to observe: 

Lemma 1.1. For any small 9 > there is an a > such that the spectrum of A is 
contained in the region 

(1.7) M e , aM = {z G C I I Imz| < 9(Rez + a), Rez> -h}. 

Proof. Let K = || Imb\\ Lcc ^y From the inequalities for a b (u, u) we see that for u G i/ 1 (0), 

I lmab(u,u)\ = I Im(&7ou,7o«)| < if (£||-u||i + C(e)||u|| ) 

< Kec^ 1 (Rea b (u,u) + &i|Mlo) + ^C( £ )IMIo 
= Kec^ 1 Rea b (u,u) + (Kec^h + KC(e))\\u\\l. 
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This (together with (1.4)) shows that for « ^ 0, a&(it, w)/||it||o has its values in M^^, 
where 6 = Kec^ can be taken arbitrarily small, a = Kec^ k\ + KC(e). The numerical 
ranges of A and A* are contained in this set, which then also contains the spectra. (More 
details for this kind of argument can be found in [G09], Sect. 12.4.) □ 

The Neumann-type boundary operator (0.2) enters in the "halfways Green's formula" 

(1.8) (Au,v) -a(u,v) = (w,7o^)l 2 (e), 

for smooth u and v. It is known e.g. from [LM68] that 71 and v extend to continuous 
mappings from if 1 (0) n D(A max ) to #~i(E), such that for u G if 1 (0) H L>(A max ), v G 
ii/" 1 (0), (1.8) holds with the scalar product over E replaced by the sesquilinear duality 
between H~^{T) and #*(£). Then 

(1.9) (Au,v)-a b (u,v) = (i/m,7 v) w _i (e) ^i (e) - (&7o«, 7o«)l 2 (e), 
and hence 

D(A) = {«£ if 1 (ft) n £>(A max ) I z/w = &70-U in fl" _ 3(E)}, 

representing the Robin condition z/tt — Iryou = 0. 

For 6 = 0, the condition is vu = 0, defining what we call the Neumann realization 
A u ; it is selfadjoint with D(A U ) C H 2 {Vt). It is well-known that when 6 is smooth, then 
D(A) C # 2 (E). 

Lemma 1.2. When b G L 00 (S), t/ie domain of A satisfies 

D(A) C Hi(Q)nD(A 

max ^ • 

Proof. When u G -D(A), then u G i/ x (0) implies 7 w G iPs(E) C P 2 (E). Multiplication 
by 6 is continuous on L 2 (E), so 670^ G L 2 (E). Then also vu = b^ou is in L 2 (E). By the 

3 — 

ellipticity of the Neumann problem, Au G L 2 (0) with vu G P 2 (E) imply u G if 2(0). □ 

When 6 has some smoothness or piecewise smoothness, we can get more regularity: It is 
known that when b is in the Bessel potential space ip^E) with r > (n — 1) /p, p > 2, then 
multiplication by 6 is continuous in iP(E) for |s| < r (cf. e.g. Johnsen [J95]). In relation 
to Holder spaces C r and Besov spaces B r there are inclusions 

(1.10) C-+ M (£) ^ B£*(E) ^ S;+/(E) ^ flJ(E), any 5 > 0, 

so also functions in these spaces preserve iiP(E) for \s\ < r. (A summary of the relevant 
facts on function spaces is given e.g. in [AGW10], Sect. 2.) Note that any e > can be 
included as an r + 26 by taking r G ]0, e[ , 5 = (e — r)/2 and p > (n — l)/r. 

When X(E) is a function space over E, we say that b is piecewise in X, when the 
(n — 1) -dimensional manifold E is a union Ei U • • • U E j of smooth subsets Ej with disjoint 
interiors (such that the interfaces are smooth (n — 2)-dimensional manifolds), and there 
are functions bj G X(E), such that b equals bj on the interior E°, for j = 1, . . . , J. 

It is well-known that multiplication by Is. is continuous on iP(E) for all \s\ < \. 
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Proposition 1.3. 

1° Let b G Hp(E) with r > (n — l)/p, p > 2 (it holds if b is in one of the spaces in 
(1.10);. Then D(A) C #i +r (fi) ifr<\, D(A) C # 2 (fi) t/r > §. 

2° Let 6 6e piecewise in H r p {Y) with r > (n - l)/p, p > 2. Then D(A) C #i +r (fi) t/ 
r < \, D(A) C H 2 ~ £ (n) for anye>Qifr>\. 

Proof. As already noted, it G implies 7o"U G -ff^(E). In the case 1°, multiplication 

by 6 preserves .£P(E) for |s| < r, so 670 it G if mm ^ r '2 ^(E). Then also vu = It/qU is in 
# min {^}(£), and now Aw G L 2 (0) with z/w G # min { r '5>(E) imply u G H'i +r (Vt) if r < J, 
u G H 2 (Q) if r > |, by the ellipticity of the Neumann problem. 

In the case 2°, since b = J2j=i^j^j: multiplication by b maps H r (E) into itself if 
r < |, and into H^~ £ , any £ > 0, if r > |. Completing the proof as under 1°, we find 
that u G H? +r (fl) if r < §, u G H 2 ~ e {9) if r > ±. □ 

Let us regard A from the point of view of the general extension theory of [G68], as 
recalled in [BGW09], [G08], [GlOa]. 

We take the Dirichlet realization A 1 as the reference operator, assumed to have a positive 
lower bound. (Seen from the point of view of [G68], [BL07] uses instead the Neumann 
realization A v as the reference operator.) The operator A corresponds, by the general 
theory of [G68], to a closed densely defined operator T: V — > W, where V and W are 
closed subsets of Z = ker A max , and D(T) is dense in V; and this in turn is carried over 
by use of the homeomorphism 70 : Z ^> H~ 2(E), to a closed operator L : X — > Y*, with 
domain D(L) dense in X, where X and Y are closed subspaces of H~ 2 (E). Here X = 70^, 
Y = 70W and D(L) = ^D{T) = <y D(A). 

Proposition 1.4. The operator L : X — » Y corresponding to A by [G68] /10s X = Y = 
H~z (E), and acte H/ce 6 — P°^ a domain contained in if 1 (E). When b is real, L is 
selfadjoint as an unbounded operator from H~^(T,) to ifi(E). 

Proof. Besides the description referred to above, we shall use the observations on operators 
defined by sesquilinear forms worked out in [G70] (and partly recalled in [G09], Ch. 13.2, 
see in particular Th. 13.19). Since the domain of a,b(u, v) equals i/ 1 (0), T is defined from 
a sesquilinear form t(z, w) with domain i7 1 (0) n Z dense in Z, and hence V = W = Z. It 
follows that X = Y = H~2(E), and L is densely defined and closed as an operator from 
H~z (E) to H 5 (E). The adjoint L* is of the same type and corresponds to A*. When b is 
real, A is selfadjoint as noted above; then L is selfadjoint. 

In the interpretation of the extension theory, A represents the boundary condition 

7 mG-D(L), Tu = L'Jqu; 

where Vu = vu — Pi? 'Jqu, so L-f u = vu — Pi? 70 it when u G D(A). (P~ „ is the operator 
mapping Dirichlet boundary values to Neumann boundary values for solutions of ( A— A)it = 
0; more on this below.) Since the functions in D(A) also satisfy vu = 670^, we see that L 
acts like 

L^=(6-P 7 °>. 
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By Lemma 1.2, D(A) C #§(£), so D(L) = 7o-D(^) C ^(E). □ 

When we replace A by A — A, where A is in the resolvent set f?(Ay) of Ay, we get for 
the corresponding operator L A : 

L A acts like 6 - P* v , with D(L X ) = D(L) C #*(£). 

For A G £»(Ay) fl £>(A), there holds a Krein resolvent formula (shown in [BGW09], Th. 
3.4): 

(1.11) (A - A)" 1 = (Ay - A)" 1 + Kfo*)- 1 ^)*- 

Here K x is the Poisson operator for the Dirichlet problem, i.e. the solution operator 
K x : (p i — y u for the problem 

(A — X)u = on O, 7ott = <p on E; 

it maps i7 s ~^(E) — >■ H S (Q) continuously for all s, and the adjoint (if A )* maps e.g. L 2 (0) 

to #3(E). 

We can use this to show a spectral estimate for (A — A) -1 — (A v — A) -1 , going via 
differences with the Dirichlet resolvent. The argumentation is not the same as that of 
[BLLLP10], which uses a Krein formula based on the Poisson operator for the Neumann 
problem. 

The spectrum of Ay is contained in a positive halfline [co, oo[ , and the spectrum of A v 
is contained in a larger halfline ] — k, oof, cf. (1.2). For A G C\] — k, oo[, the Dirichlet- 
to-Neumann operator P x = vK x is a homeomorphism from ii" s (E) to ii" s_1 (E) for all 
s6i, with inverse P x tl , the Neumann-to-Dirichlet operator. Then we can write 

(1-12) L\ = (b- P x „)<p = (bP*„ - l)P* v <p, for <p G D(L). 

Since P x is of order —1, it is compact in L 2 (E). Then 6P A 7 — 1 is a Fredholm operator 
in L 2 (E), as noted also in [BLLLP10]. If A is such that: (1) L x is invertible (from D(L) 
to (E)), (2) bP^ - 1 is invertible in L 2 (E), then the inverse of L x must coincide with 
the inverse of (bP^ - l)P x v on H^{T). 

For bP x n - 1, we get invertibility as follows: We have as a simple application of the 
principles in [G96] (cf. Th. 2.5.6, (A.25-26)) that 

ll-P^iz^H-ff^CE) — IMIh s +i,m(e)> ||y||fl-»- 1 ."(s) — II^^H/z^cs), 

uniformly in = | A j 2 for A — > 00 on rays in C \ R+; this holds since P x v is parameter- 
elliptic of order 1 and regularity +00 on the rays in C\R+. In particular, one has on such 
a ray {A = fi 2 e lv } with 77 G ]0, 2n[ , for s G [0, 1] and \i > 1, 

II-P^IIh-Ce) + (^>1^ A 7 </?IIl 2 (e) < Cmin{||^|| Hs -i (E) , (^> s_1 |MIl 2 (e)}, 
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so the norm of P A 7 in L 2 (E) is O^) -1 ) on the ray. Take /iq so large that \\bP x 1 \\c(L 2 (T:)) < 
5 < 1 for fi > (Uo, then bP x n — 1 is invertible as an operator in L2(E) for fi > ^o, with a 
bounded inverse ipP^ — l) -1 : 

oo 

(1.13) (6P* 7 - l)" 1 = -1 - 5>P„, 7 ) fc , converging in £(L 2 (E)). 

fc=i 

Then 6 - P^ has an inverse 

(1-14) (b - P 7 A , J" 1 = P,, 7 (W^ 7 - l)" 1 . 

For L x we know from the extension theory that L x is bijective from D(L) to P"^(E) if 
and only if A G g(A). It follows from Lemma 1.1 by a simple geometric consideration that 
for each ray {A = fj, 2 e ir1 } with rf g]0,27t[, there is a fi\ such that such that A G g(A) for 

fl > fix. 

For /i > maxj^o, /Ui}, both (1) and (2) are satisfied, so then 
(1-15) (L")- 1 = (b- P^)- 1 = P^(bP^ - l)" 1 on PT* (E). 

We note in particular that 

(1.16) D(L X ) = & G Pf^E) | (6 - P 7 A 5 > G *(E)}, 

for such A. Now D(L) = D(L X ), and - P^ is bounded from P"^(E) to #*(E) (cf. 
[BGW09], Rem. 3.2), so we conclude that 

(1.17) £>(L) = G P^E) | (6 - P 7 °> G ff*(E)}. 

It follows moreover that (1.16) holds for all A G £>(Ay). 
This shows the main part of: 

Theorem 1.5. The domain of L satisfies (1.17), and it is also described by (1.16) for any 
A G q{A 1 ). 

On each ray in C\1R+, A is in g(A) and (1.15) holds for |A| sufficiently large. For such 

(1-18) (A - A)" 1 - (A, - A)" 1 = K x P u , 7 (bP^ - l)- 1 ^)* . 

Proof. The statements before formula (1.18) were accounted for above, and the formula 
follows by insertion of (1.15) in (1.11). □ 
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2. Spectral estimates. 

Spectral estimates for resolvent differences will now be studied. A classical reference for 
the basic concepts is the book of Gohberg and Krein [GK69]; some particularly relevant 
facts were collected in [G84], supplied with additional results. We shall include a short 
summary here: 

For p > 0, the space C p is the Schatten class of compact linear operators T (in a Hilbert 
space H) with singular value sequences (sj(T)) j€ N £ and & p denotes the quasi-normed 
space of compact operators T with Sj(T) = here & p C C p+£ for all e > 0. 

The rules shown by Ky Fan [F51] 

(2.1) Sj+k^iT + T') < Sj (T) + s k (T'), Sj+fc _ 1 (TT / ) < s i (T)s fc (T / ), 
imply that C p and & p are vector spaces, and that a product rule holds: 

(2.2) 6 p • & q C 6l/(p-i+g-i), C p • C g C Ci/( p -i +(? -i). 

Moreover, the rule 

(2.3) fli (ATB) < P|| S ,-(T)||B|| 

implies that © p and C p are preserved under compositions with bounded operators. We 
mention two perturbation results: 

Lemma 2.1. 

1° // Sj(T)j 1 / p -> C and s^T')] 1 ^ ->■ /or j ->■ oo, tfien Sj -(T + T')j 1/p ->■ C /or 
J ->■ oo. 

2° If T = T M + T' M for each M G N, w/iere Sj{T M )j 1/p ->■ C M /or j ->■ oo and 
Sj(T' M )j 1/p < e M for j G N, witfc C M ->■ Co and e M ^ for M ^ oo, then s i (T)j 1 / p ->• 
C /or j ->■ oo. 

The statement in 1° is the Weyl-Ky Fan theorem (cf. e.g. [GK69] Th. II 2.3), and 2° is 
a refinement shown in [G84], Lemma 4.2.2°. 

We also recall that when S is a compact n'-dimensional smooth manifold (possibly with 
boundary) and T is a bounded linear operator from L 2 (£) to f/"*(H) for some £ > 0, then 
T G @ n '/t as an operator in L 2 (S), with 

(2-4) sjiT)^' <C\\T\\ ciL2 , Ht) , 

C depending only on S and t. See [G84], Lemma 4.4ff. for references. 

The Poisson operator is continuous from H S ~^(T,) to H S (Q) for all s G R, and its 
adjoint is a trace operator of class and order —1 in the pseudodifferential boundary 
operator calculus, hence is continuous from H S (Q) to H S+ ^(Y,) for s > —\. Then the 
composition K^*K^ is continuous from L 2 (E) to if 1 (E), so in view of (2.4), K^*K^ G 
& n -i and hence G @( n -i)/(i/2)) as operators in L 2 (E). The singular numbers of K** 
have the same behavior. Moreover, since is a pseudodifferential operator of order — 1 
on S, it lies in (E5 n _i when considered as an operator in L 2 (E). 
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Theorem 2.2. Let 6 G L OC3 (E). For any A G o(A) n q(A v ), 
(2.5) (1 - A)" 1 - (A v - A)" 1 G e (n _ 1)/3 . 

Proof. First assume that A lies so far out on a ray in C\R+ that the statements in Theorem 
1.5 are valid. 

Applying (1.18) to our A and also to the case 6 = (the Neumann realization), we find 
by subtraction: 

(A - A)" 1 - (A u - A) -1 = (A - A)" 1 - (A, - A)" 1 - ((A, - A)" 1 - (A 1 - A)" 1 ) 
(2-6) = K^[(bP^ - l)" 1 + 1]K? 

= K$P*„(bP?„ - l)~ l bP^Kf. 

The last expression is composed of the operator in 6( n -i)/(i/2)) the adjoint of with 
the same property, two factors in <5 n _i and the bounded operators (bP£ n — 1) _1 and 
6, so it belongs to 6( n -i)/3> by (2.2). 

Now let A' be an arbitrary number in g(A) n o(^4^). We use the following refined 
resolvent identity as in [BLLLP10]: 

(2.7) (S - A') -1 — (T — A') -1 

= (1 + (A' - A)(T - A')- 1 )((5 - A)" 1 — (T — A)- X )(l + (A' - A)(5 - A')" 1 ), 

valid for A, A' G o(T) fl g(S). Applying it to S = A and T = A v for A as above and 
A' G g(A) fl £(A„), we find that (A — A') -1 — (A^ — A') -1 is a composition of an operator 
in <5( n _i)/3 with two bounded operators, hence lies in @( n _i)/3, as was to be shown. □ 

There is an obvious corollary: 

Corollary 2.3. Let 61,62 G L 00 (E), and denote the corresponding realizations of Robin 
conditions vu = 6170W resp. vu = 6270^ by A\ resp. A^. For any A G g{A\) fl 0(^2), 

(2.8) (A 1 - A)" 1 - (1 2 - A)" 1 G e (n _ 1)/3 . 

Proo/. Write (Ax - A) -1 - (A 2 - A) -1 as the difference between (Ax - A) -1 - (A v - A) -1 
and (A 2 — A) -1 — (A 7 — A) -1 , then the result follows from Theorem 2.2 (and (2.7)) since 
& p is a vector space. □ 

Formula (1.18) also allows us to show a spectral asymptotics estimate for (A — A) -1 — 
(A 1 - A)" 1 that was obtained in the smooth case for selfadjoint realizations and negative 
A in Grubb [G74], Sect. 8, and Birman and Solomiak [BS80]. In the former paper it is 
shown, also for 2m-order problems, that the operator is, on the orthogonal complement of 
its nullspace, isometric to an elliptic pseudodifferential operator on E of order —2m (which 
has the asserted spectral asymptotics); in the latter paper exterior domains are included. 
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Theorem 2.4. Let b E L OC3 (E). For any A G o(A) n g(A^), 

(2.9) - A)" 1 - (A, - A)" V /( " _1) C 2/(n - 1} /or j oo, 
where Co is the same constant as in the case 6 = (where A = A v ), namely 

(2.10) Co = (w _ 1)( L )n -i / / _ (UPlU^M^bT 72 )"- 1 ^^)^'; 

here k°(x' , x n , £') is £/ie principal symbol-kernel of andp°(x', £') is the principal symbol 
of P x 

Proof. Since the details are perhaps not very well known, we first give a proof of (2.9)- 

(2.10) in the case 6 = 0. We have as an easy special case of (1.18) that 

(2.11) (A u - A)" 1 - {A 1 - A)" 1 = -K*P v „K>f = G v . 

This is a singular Green operator with principal boundary symbol operator 

g°Ax',e,D n ) = -k o (x>,e,D n )p°(x',e)k o (x',e,D n y 

in local coordinates, where k° and p° are the (A-independent) principal symbols of K x 
and P* r At each k°(x' , D n ):C ->■ L 2 (R+) maps v E C to fc°(x', x n , £' )v, 

where k°(x', x n , £') E <S(R + ) is the symbol-kernel. In the case A = —A it equals e - ^ x ™, 
and it has a similar structure for general A (cf. e.g. [GS01], Sect. 2.d). The operator 
k°(x', £>n) *: ^ 2 (K+) -)• C maps u(x n ) to (u, A;°)l 2 (h + ). Thus k°(x', £>n)*/c o 0r', £' , A») 
is the multiplication by ||&°||| 2 ( R+ )) and k°(x\ D n )k°(x', -D n )* is the rank 1 operator 
mapping it to (u,k°)k°. The latter operator has the sole eigenvector /cj 1 = fe /||fc°|| with 
a positive eigenvalue ||/c°|| 2 (besides eigenvectors in the nullspace), so its trace equals the 
eigenvalue. The middle factor p° is just multiplication by a scalar; for A = —A, it equals 

-irr 1 - 

By [G84], Th. 4.10, since G v is a singular Green operator of order —2 and class 0, 

(2.12) Sj (G u )j 2 ^-V _> C(<£) 2 /<»-i) for j oo, 
where 

(2.13) C(gg) = (w _ 1)( L)~-i / / tr^a^x',^^)^^,^^))^" 1 ^ 4 ]^^)^. 
Here 

2°(^^D n )*2°(a/, £',£>„) 

= ||£ V, x B , O ||| 2(K+) bV, 01 2 *V. AO* V. AO*- 
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This is a rank 1 operator with eigenvalue \\k ||£ 2 |p | , so 

trM(x',e,D n r g ^x',^,D n )) {n - 1)/4 ] = (fk°\\ L2 \p°\ 1/2 ) n -\ 
and (2.10) follows. 

Now the case of general 6: For large A on rays in C \ R+ as in Theorem 1.5 we write 
formula (1.13) as 



k=0 



(2.14) (bP^ - l)" 1 = -1 - bP^S, where S = ^(6P, A 7 ) fc G C(L 2 (E)). 

Then we have from (1.18): 
(2.15) 



(A - A)" 1 - (A 7 - A)" 1 = K*P v „{-\ - bP^S)K* 



The first term equals (A U — X)~ 1 — (A 1 — X)~ 1 and satisfies the spectral asymptotics estimate 
(2.9) with (2.10). The second term is in <5( n _i)/3, in view of the mapping properties of its 
factors, as in the proof of Theorem 2.2. By Lemma 2.1.1°, it follows that the sum of the 
two terms has the asymptotic behavior (2.9). 

General A G g(A) n g(A^) are included by use of the resolvent identity (2.7), which 
gives the operator as a sum of a term with the behavior (2.9) and terms in ©( n _ 1 - ) /( 2 +t) 
with t > 0, using that (A 7 - A) -1 G 6 n / 2 and (A - A) -1 G 6 n /(3/2)- Then Lemma 2.1.1° 
applies to show (2.9) for the sum. □ 

Spectral asymptotics estimates for the resolvent difference (2.5) are harder to get at, 
since b here enters in the principal part of the operator. However, with a little smoothness 
of b we can obtain the spectral estimate by reduction to a case that allows an approximation 
procedure. 

We consider the resolvent difference of two general Robin problems from the start, since 
the asymptotic property is not in general additive. 

Theorem 2.5. Assume that 61,62 G ifp(E), where r > and p > (n — l)/r, p > 2; this 
holds if the bi are in one of the spaces in (1.10). Define Ai as in Corollary 2.3. Then for 
A G g(A 1 )ng(A 2 ), 

(2.16) BjdM - A)" 1 - (A 2 - A)- 1 )//^- 1 ) -> C(<7 ) 3/(n_1) for j -+ 00, 
where 



(2-!7) C(g°) = / / (||P||! 2(R+) b°| 2 |6 2 -6 1 |)^- 1 )/ 3 rfa;(n^. 



(n-l)(27r)«-i y E> / !e| = 1 VI1 HMR+) 

Proof. First let A be large on a ray in C \ M_|_ such that Theorem 1.5 applies to A\ and 
A 2 . Using (2.14) in the form 

(biP}„ - I)" 1 = -1 - b t P}„ - {hP^fSi 
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we have that 

(hP^ - l)" 1 - (b 2 P^ - l)" 1 = (6a - h)P^ - (b 1 P^) 2 S 1 - (b 2 P^) 2 S 2 . 

Then we get using (2.6): 
( 2 - 18 i 

( Al - A)" 1 - (A 2 - A)" 1 = (A 1 - A)" 1 - (A v - A)" 1 - ((A 2 - A)" 1 - (A v - A)" 1 ) 

= K^P^[(hP^ - l)" 1 + l]K~f - K*P^[(b 2 P^ - l)" 1 + l]Kf 

= K*P^{b 2 - h)Pi„K? - K$P u x 1 (b 1 P^) 2 S 1 Kf + K^P^{b 2 P^fS 2 K^ 
= G + F 1 + F 2 . 

In the terms Fi we use for one of the factors hP^ that bi preserves iP(E) for \s\ < r (see 
the text before Proposition 1.3), so that &i-P^ 7 maps L 2 (E) continuously into H r ' r' = 
min{r, 1}. So this factor is in <5( n _i)/ r ', together with the usual two factors in @( n -i)/(i/2) 
and two factors in @ n _i, whereby the full composed operator Fj is in (5( n _i)/( 3+r /). It 
will not influence the spectral asymptotics. 

In the term G, let us denote b 2 — b\ = b. We write b for each M G N as a sum 

(2.19) b = b M + b' M , 

where b M G C°°(E) and sup x /g^ |6^(x')| < 1/M; this is possible since 6 is continuous on 
the smooth compact manifold E. Accordingly, we write G = Gm + G' M with 

Gm = KyP^jbMP^-yK^ , G' M = K 1 P vl b' M P vl K 1 . 

Here G' M is a composition of fixed operators with the usual (5 p -properties and a factor 
b' M whose norm in £(L 2 (E)) is < 1/M; this implies that 

(2.20) svLp Sj (G' M )ftt n -V < C/M, all M, 

3 

for a suitable constant C, in view of (2.3). 

The term Gm is treated by application of the tools in [G84]. Since 6m G C°°, Gm 
is a genuine singular Green operator of order —3 and class 0, with polyhomogeneous 
symbol. The principal symbol g M is the symbol of the boundary symbol operator (in local 
coordinates) 

(2.21) g° M (x', D n ) = k°(x', D n )p°(x', OM*>V. £')*V. £>„)*. 
It follows from [G84], Th. 4.10, that 

(2.22) Sj (G M )f/^ CG&) 3/(b_1) for j oo, 

where 
(2.23) 

C(9 M ) = (n-D^n-i / / trf^^,^,^)*^^,^,^))^-^ 6 ]^^)^. 
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As in the analysis of g„*g® in the proof of Theorem 2.4, now with the middle factor p° 
replaced by p°bMP°, we find that 

tr [( 9 Ux',e,D n y g Ux',t;',D n )) {n - 1)/6 ] = (||^i>TlM 2 ) (n - 1)/6 

= (ll^ ll! 2 b°l>M|) ( - 1)/3 , 

and hence 



(2.24) C(g° M ) = (w _ 1)( L)n-i ^II^IIL^)^! 2 ^^!)^- 1 ^ 3 ^^^- 

When M — > oo, 6m(^') — > b(x') uniformly in x' , so 

(2.25) C(g° M )^C(g°), where 

c(g°) = (w _ 1)( L)~-i j[/ #=i (ll*°lli a( R,)b l 2 |6|) (n - 1)/3 ^(O^ 



with b = b 2 — b\. 



Now we first apply Lemma 2.1.2° to the decompositions G = Cm + this shows 
that G has the spectral behavior in (2.16). When F\ and F 2 are added to G, we can use 
Lemma 2.1.1° to conclude that also G + Pi + F 2 has the spectral behavior in (2.16). 

Finally, general A G g(A) fl g{A v ) are included by use of the resolvent formula (2.7) as 
in the preceding proof. □ 

In the case A = —A, where /c° and p° are independent of x', the formula for C(g°) 
reduces to a constant times J s \b 2 — b\\ ( ' n ~ 1 ^' i dx' . 

3. Coefficients with jumps. 

It possible to extend the result of Theorem 2.5 to cases where b has jump discontinuities, 
by use of special results for pseudodifferential operators (from here on abbreviated to 
V'do's). In showing this, we also supply the general knowledge on spectral asymptotics for 
V'do's multiplied with nonsmooth functions. 

Let S be a compact n'-dimensional C°°-manifold without boundary, and assume that 
it is divided by a smooth (n' — l)-dimensional hypersurface into two subsets and S_ 
(n'-dimensional C°°-manifolds with boundary) such that S = U H_, S ^ fl S°_ = 0, 
9S + = <9£_. (Since the sets need not be connected, this covers the situation of J smooth 
subsets described before Proposition 1.3.) We denote by r ± the restrictions from H to H±, 
and by e ± the extension-by-zero operators from functions on S± to functions on S: 

, f u on E± 

e u — ' 



on S T . 

Multiplication by the characteristic function lg for S + can also be written e + r + ; similarly 
1=_ = e~r~ . 

It is well-known (as recalled e.g. in [G84], Lemma 4.5) that when P is an iV x iV-matrix 
formed classical i/>do on S of negative order —t, then it satisfies the spectral asymptotics 
formulas for j — > oo: 

Sj(P)j t/n ' -)• C(p°) t/n ' in general, 
^'^ ±\f(P)j t/n ' -)• ^(p )*/"' if P is selfadjoint, 
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where, respectively, 

C(P°) = [ I _ tr [{p°(x,0V(x,0f /2t } MOdx, 

(3.2) 77 l=1 

Let us also recall the result of Laptev [L77, L81]: 

Proposition 3.1. Let P be a classical pseudodifferential operator on S of negative order 
-t. Then G 6 K -i )/t . 

(Expressed in local coordinates, this means that the operator whose kernel is the re- 
striction of the kernel of P to the second or fourth quadrant, picks up the boundary 
dimension in its spectral behavior. For -0do's having the transmission property at <9£ + , 
this is confirmed by the results of [G84].) 

The rules in the following are valid also for N x iV-matrix formed operators P and 
factors 6, and would then need a trace indication tr in the integrals; we leave this aspect 
out here for simplicity. 

Theorem 3.2. Let P be a classical pseudodifferential operator of negative order —t, such 
that (Pu,u) > for u G 1/2(5). Then P(+) = 1e: + P1e: + satisfies the spectral asymptotics 
formula 

(3.3) «i(P(+))j t/B ' <P(+)) t/n ' for J -X oo, 
where 



(3.4) 



< p w) = ^hw I j^_y{x^yp\x^)) n ' /2t du{i)dx 

[ [ p°(x,O n,/t MOdx. 
Je, J\e\=i 



ri(2ir) r 



Proof. The principal symbol p° is > 0; which explains the second identity in (3.4). Intro- 
duce two C°° cutoff functions Q\ and (2 taking values in [0, 1] such that Q\ = 1 on and 
vanishes outside a neighborhood of H+, and (2 = on £_ and is 1 outside a neighborhood 
of H_. We shall then compare P( + ) with the operators (all are compact in 1*2(5)) 

Pi = C1PC1 and P 2 = C2PC2. 
When u G £2(2), denote e ± r ± u = u±. We have for Pi, since (iu + = u + : 
(Piu,u) = (Piu+,u+) + (PiU+,U_) + (PiU_,U+) + (Piti_,U_) 

= (P(+)«> w ) + u ) + (PCiw-> &«-)» 

where P = 1h_PiPs + + ls + Piis_. Since Pi is a classical -i/'do of order —t on S, it has the 
spectral behavior in (3.1)-(3.2) with the limit C(pi)*/ n '; here 

C ^) = ^itw I I (bp°(x,t)Ci) n,/t Mt)<k- 
n[27T > y SU p P ci^i?i=i 
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Moreover, P is of the type considered in Proposition 3.1, hence lies in &(n>-i)/t- Then 
by Lemma 2.1.1°, Pi — R likewise has the spectral behavior in (3.1)-(3.2) with the limit 
C(piY^ n . Now observe that since P is nonnegative, (P(iU-,(i u -) > for all u G L 2 (r,). 
Thus we have: 

(3.5) (P(+)U, u) < ((Pi - R)u, u), for all u G L 2 (E). 

Both operators P(+) and Pi — P are selfadjoint nonnegative, so the s-numbers are the same 
as the eigenvalues, and the minimum-maximum principle implies in view of (3.5) that 

(3.6) Sj(P (+) ) < - R), for all j. 
It then follows from the limit property of the Sj (Pi — R) that 

(3.7) limsup^^lV))//"' < C(pS)*/»'. 
For the comparison with P2 we write, using that ( 2 u+ = ( 2 u, 

(P {+) U, U) = (( 2 P( 2 U+, U+) + ((1 - C2)P(1 " C2)«+, U+) + ((1 - C2)i 3 C2«+, «+) 

+ (C 2J P(i-C2)«+,«+) 

> (( 2 P( 2 U, U) + ((1 - C 2 )PC2U, U+) + (C 2 P(1 " C 2 )«+, «) 

= ((C2PC2 + (1 - C 2 )PC 2 + C 2 P(i - C 2 )K m) + «), 

where Pi is a sum of terms as in Proposition 3.1. Then since Sj(P^) = Xj(P^) > 

A+(C 2 PC 2 + (i - C 2 )PC 2 + C 2 P(i - C2) + Pi), 

(3.8) hminf^ ooS ,(P (+) )//n' > C+(( 2 p°( 2 + (1 - ( 2 )p°( 2 + ( 2 p°(l - C 2 )Y /n ' • 

Since C(p\) and C + (Q 2 p°Q 2 + (l — Q 2 )p°Q 2 + Q 2 p {l — ( l2 )) come arbitrarily close to c(P( + )) 
when the support of Ci shrinks towards S + and the support of 1 — £2 shrinks towards H_, 
we conclude that (3.3) with (3.4) holds. □ 

This leads to a result on compositions of ^do's with discontinuous factors, which seems 
to have an interest in itself: 

Theorem 3.3. Let P be an operator composed of I classical pseudodifferential operators 
Pi, . . . , Pi of negative orders —ti, . . . , — ti and l + l functions bi, ... , bi+i that are piecewise 
continuous on S with possible jumps at <9£_|_ (so the bk extend to continous funcxtions on 
S_|_ and on z,-); 

(3.9) P = b 1 P 1 ...hPib l+1 . 



Let t = ti + 
(3.10) 



• + ti . Then P has the spectral behavior: 

SjiP)^' -+ c(Py/ n ' for j -> 00, 
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where 



(3-11) • b 1 (x)p° 1 (x, . . OW*))" 7 * 

|6l.. 6 I+1 p?...p?r , /t da ,^) da; . 



n'(2-n-y 




l€l=i 



Proof. We can write 



P*P = b l+1 P l * ...P*b 1 b 1 P 1 ...P l b l = l~ + P*Pls + +1~_P*P1~_ +R, 

where i? = l~ + &/+i-Pf • • • P/M~_ + l~_fy+i-Pf • • • -P/M-+- Inserting 1 = l s+ + l s _ at 
each factor b\~ or bk in R and multiplying out, we obtain it as a sum of terms of order —t, 
each containing at least one factor of the type in Proposition 3.1. Thus R G & n > /(t+8) with 
a 5 > 0. For the term 1= + P* Pl~. + , we proceed as in Theorem 2.5. We can assume that 
bk is extended from S + to a continuous function bk on S. Each bk is approximated by a 
uniformly convergent sequence bkM of C°°-functions on S. For each M, 

PmPm = &/+i,m-P/* • • • P\b\M.b\MP\ ■ ■ ■ blMPlbl+l,M 

is a classical nonnegative V'do of order — t, so Theorem 3.2 applies to the operator with 
1= + before and after, and gives the corresponding spectral asymptotics formula. Since 
PmPm — P*P can be written as a sum of terms where each has a small factor bkM — bk 
or bkM — bk, we have for M — > oo that 

(3.12) sups J -(l s+ PjE f P M ls + - l~ + P*Pl~ + )j t/n ' -)• 0. 

Then Lemma 2.1.2° implies a spectral asymptotics formula for 1h + -P*-P1e: + , with the 
constant as in (3.11) but integrated over H + . — There is a similar result for 1s_-P*-P1e:_ , 
relative to H_. 

Now since Z/2(H) identifies with the orthogonal sum of L2CE.+) and L2(S_), the spectra 
are simply superposed when the operators are added together. The statement Xj (T)j t ' n — > 
c{T) t l n ' for j ->■ 00 is equivalent with N'(a;T)a n ' /l ->■ c(T) for a ->■ 00, where N'(a;T) 
counts the number of eigenvalues in [1/a, 00 [ ; superposition of the spectra means addition 
of the counting functions. (More on counting functions e.g. in [G96], Sect. A. 6.) Thus 
1= + P*P1=: + + 1s_P*P1e:_ has a spectral asymptotics behavior where the constant is 
obtained by adding the integrals for l~. + P*P^-~. + and 1s_P*P1e:_, so it is as described 
in (3.9)-(3.11). By Lemma 2.1.1°, the behavior keeps this form when we add R to the 
operator. □ 

A similar theorem holds for matrix formed operators Pk and factors with c(P) 
defined by the first expression in (3.11); here of course it cannot be reduced to the second 
expression unless all the factors commute. 

A special case of the situation in Theorem 3.3 is the case of bP, where P is a classical 
V'do and b is a piecewise continuous function. We need a case with interspersed factors bk 
in our application below. 

We can now show: 
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Theorem 3.4. The conclusion of Theorem 2.5 holds also when b\ and b 2 ore piecewise in 
Hp(E) for some r > as in Theorem 2.5, b 2 ~ &i having jumps at a smooth hypersurface. 

Proof. We use again the decomposition in (2.18): 

(Ii - A)" 1 - (A 2 - A)" 1 = G + F 1 + F 2 , with G = K*P^(b 2 - b 1 )P^K$\ 
F 1 = -K^P^(b 1 P u x 1 ) 2 S 1 K^\ F 2 = K^P^(b 2 P^) 2 S 2 K x \ 

and F\ and F 2 are handled as after (2.18), using that biP* maps L 2 (E) into H r (E), 
r' = min{r, | — e}. Then they are in (5( n _i)/( 3+r /). We denote again b 2 — bi = b. 

For G we proceed as follows: Let A be large negative, so that Theorem 1.5 holds. Since 
A is real, = and is selfadjoint. The j-th eigenvalue of G*G satisfies 

Xj(G*G) = XAK^bP^K^K^bP^K^ ). 

Here K x * K x equals a selfadjoint tfjdo P± of order —1; it is nonnegative on L 2 (E) and 
injective, since K x is injective: 



(^,¥>)l 2 (e) = (K^K^tp)^ = ||^||| 2(n) > c\\ip\\ 2 H _ h 



(E)' 



hence elliptic. It follows from Seeley [S67] that P\ has a squareroot P 2 = P± which is a 
classical elliptic ^do of order — |. Then we find, applying the general formula 

(3.13) X j {TT') = X j {T'T), 

with T = K$P^bP^P 2 , T' = P 2 P*„bP^Kf, that 

Xj(G*G) = X^P^bP^P^P^bP^K^) 
= X^P^bP^K^K^P^bP^P,) 
= XjiP.P^bP^P^bP^Pi). 

The operator Q = P 2 P x rf bP x 1 P\P x 1 bP x ri P 2 is an operator to which Theorem 3.4 applies, 
and it gives a spectral asymptotics formula with the constant defined as in (3.11), with 
n' = n — 1. Since pi = , p 2 = ||&°||l 2 5 the formula can be rewritten in the form 

(2.17). 

The proof is now completed in the same way as in the proof of Theorem 2.5. □ 

The results can be extended to exterior domains by the method of [G10]. 

In a forthcoming paper we shall treat the question of spectral asymptotics for the mixed 
problem for —A + ao, where the boundary condition jumps from a Dirichlet condition to 
a Neumann condition at a smooth hypersurface of E. Here we moreover need to draw on 
the analyses of nonstandard pseudodifferential operators, as in Shamir [S63], Eskin [E81], 
Birman and Solomiak [BS77] and many later works. 
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